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Neutrino processes in semi-transparent supernova matter, opaque to semi- 
transparent protoneutron star matter, and catalyzed neutron stars are discussed. 



Q ■ 1. Introduction 

H ■ 

' In recent years, the study of neutrino emission, scattering, and absorption in 

. matter at high density and/or temperature has gained prominence largely 

due to its importance in a wide range of astrophysical phenomena. Energy 
L J ■ loss in degenerate helium cores of red giant stars, ^'^ cooling in pre-white 

. dwarf interiors,'^ the short- and long-term cooling of neutron stars, ^'^ the 

' deflagration stages of white dwarfs which may lead to type la supernovae,^'^ 

explosive stages of type II (core-collapse) supernovae,* and thermal emis- 
sion in accretion disks of gamma-ray bursters,^' are examples in which 
neutral and charged current weak interaction processes that involve neu- 
trinos play a significant role. Since this is a vast subject, we will highlight 
some recent developments in the context of core-collapse supernovae and 
neutron stars from their birth to old age. 

2. Supernovae 

In unravelling the mechanism by which a type-II supernova explodes, the 
implementation of accurate neutrino transport has been realized to be 
critical. The basic microphysical inputs of accurate neutrino transport 
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coupled in hydrodynamical situations are the differential neutrino produc- 
tion and absorption rates and their associated emissivities. The processes 
and precise forms in which such inputs are required for multienergy treat- 
ment of neutrinos for both sub-nuclear and super-nuclear densities (nuclear 
density po ~ 2.65 x 10"^^ g cm~"^) are detailed in Rcfs. [12, 13]. 
At sub-nuclear densities, the relevant processes are: 



pair production 
the photo — neutrino process 
the plasma process 
nucleon — nucleon bremsstrahlung 
v — flavor production 



e + 7 



7 ^ t 
{n,p) 

Vi + Vi 



{n,p) + v + 1/ . 
y Vj + Vj . 



(1) 
(2) 
(3) 

and (4) 
(5) 



The relative importance of these processes depends on the temperature and 
density of ambient matter and is sketched in Figure 1. 




Figure 1. Regions of temperature and density in which the various neutrino emitting 
processes are operative. 
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Additional neutrino processes in the supernova environment include 

p + e- ^n + Ue, {A,Z) + e- ^ {A,Z -l) + Ve, (6) 

v + {A,Z)^y + {A,Z), (7) 

v + e-^v + e-, u+{A,Z) ^u+{A,Z)* , (8) 

iA,Z)* ^{A,Z) + u + D. (9) 

Reactions (6) begin the process of neutronization and decrease of lepton 
number per baryon Yi,, whose value after z^-trapping determines the masses 
of the homologous core and initial PNS, and thus the available energy for 
the shock and subsequent neutrino emissions. The equation of state also 
influences these quantities, most importantly through the nuclear symmetry 
energy. In the submiclcar density regime, the coherent scattering reaction 
(7) from nuclei in a lattice is the most important opacity source. The 
reactions (8) are important in changing the neutrino energy and in achieving 
thermodynamic equilibrium. The reactions (1) and (9) are also important in 
achieving thermodynamic equilibrium. The bremsstrahlung and modified 
Urea (n+p — > n+n + e'^ + v+v) processes dominate in many circumstances. 
For example, the production and thermalization of and r neutrinos, which 
receives contributions from reactions (1) through (5), and (9), is dominated 
by nucleon bremsstrahlung (4) for number density n > 0.005 fm~^ and 
temperature T < 15 MeV. The modified Urea process dominates the cooling 
of protoneutron stars if direct Urea processes involving nucleons, hyperons 
or other strange particles do not occur. 



3. Kernels for Neutrino Transport Calculations 

The evolution of the neutrino distribution function /, generally described 
by the Boltzmann transport equation in conjunction with hydrodynamical 
equations of motion together with baryon and lepton number conservation 

equations, is 

= BeaU) + B^Esif) + BMf) + BTp{f) . (10) 

Here, is the force acting on the particle and we have ignored general 
relativistic effects for simplicity (see, for example, Ref. [12] for full details). 
The right hand side of the above equation is the neutrino source term in 
which, BeaH) incorporates neutrino emission and absorption processes. 
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BnesH) accounts for the neutrino-electron scattering process, BuN{f) in- 
cludes scattering of neutrinos off nucleons and nuclei, and BTp{ f) considers 
the thermal production and absorption of neutrino-antineutrino pairs. 

Till recently, detailed differential information was not available for the 
plasma and photoneutrino processes. In prior works in which the total 
rates and emissivities for these processes were computed, the energy and 
angular dependences of the emitted neutrinos with 4-momenta q and q' 
were eliminated by using Lenard's identity: 

/ l^l^:'^'^^* " " '^'^'^''^"^ = ^e{qf){2qi^q'^ + q^gn ■ (H) 

Although the use of this identity simplifies considerably the calculation of 
the total emissivity, differential information about the neutrinos is entirely 
lost. On the other hand, calculations of differential rates and emissivities, 
such as 

and ^ , , (12) 



dEq dEqi d{cOS 6qqi ) dEq dEqi d{cOS 6qqi ) ' 

where Oqqi is the angle between the neutrino pairs, entail the calculation 
of the relevant squared matrix elements hitherto bypassed in obtaining the 
total rates and emissivities. Realizing this, the squared matrix elements for 
the plasma and photoneutrino processes were computed in Refs. [16, 17]. 
Some results from these recent works are highlighted below. 

3.1. The Plasma Process 

As is well known, e+e~ pairs in a plasma cause the photon to acquire an 
effective mass, which arises from electromagnetic interactions (cf. Refs. 
[16, 18] and references therein). Therefore, we can consider the photon to 
be a massive spin-1 particle that couples to the vv pair through the two 
one-loop diagrams shown in Figure 2. 

Suppressing the dependencies on {r,t) for notational simplicity, the 
source term for the plasma process can be written as 

-1 /"oo /•! /'27r 

B{f) = [1 - f]77r^ / El dE2 / d^i2 / d<l>2BF{EuE2,cos6) [1 - /] 

-I pOO pi p27T 

-fj^J^ EidE^J ^dfx^J^ d(t>2R''{Ei,E2,cosO)f, (13) 

where the first and the second terms correspond to the source (neutrino 
gain) and sink (neutrino loss) terms, respectively. Ei and E2 arc the en- 
ergies of the neutrino and antineutrino, respectively. Angular variables 
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Figure 2. Leading order Feynman diagrams describing the emission of a neutrino pair 
from the plasma process. The charged current process in which the Vl^-boson is exchanged 
produces only Vf,v^, while that in which the neutral Z-boson is exchanged results in pairs 
of all three neutrino (e, and r) flavors. 



pLi = COS 9i and 0i (z = 1,2) are defined with respect to the z-axis that is 
locally set parallel to the outgoing radial vector r. The angle 9 between the 
neutrino and antineutrino pair is related to 6i and 02 through 



cos 61 = HiH2 + Y^(l - /Uf)(l - nl) COS(01 - (p2) ■ 



(14) 



Notice that / = f{Ei,^i) and / = /(£^2,M2)- The production and absorp- 
tion kernels are given by 



Rl{E,,E2,cose) = J 



-(fc) 



SijEl E2 



S\K-Q^-Q2){2n)'{\Mn, (15) 



where the subscript Y stands for T- "transverse" or //-"longitudinal". The 
factor ^ accounts for the spin avaraging; ^ = 2 for the transverse, axial and 
mixed cases, while for the longidudinal case ^ = 1. 

The angular dependences in the kernels ii» {Ei,E2,cos 9) are often ex- 
pressed in terms of Legendre polynomials as 



R"^{E^,E2,cos9)=Y, ^} {El ,E2)Pi (cos 0) ., 



(16) 



1=0 



where the Legendre coefficients {Ei,E2) depend exclusively on energies. 

Prom Eq. (15), it is evident that the kernels are related to the neu- 
trino rates and emissivities. We first consider the production kernel 
R^{Ei,E2,cos9). The corresponding analysis for the absorption kernel 
R"{Ei,E2, cos 6) can be made along the same lines, but with the difference 
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that ub is replaced by I + ub- The neutrino production rate is given by 



which defines the kernel R^{Ei,E2, cos 9) and is to be identified with that 
in Eq. (15). The emissivity Q can also be cast in terms of BP using 



The Lcgcndrc coefficients "f"^ for the transverse part of the production 
kernels are shown in Figure 3 for I = through 3. The longitudinal compo- 
nent becomes comparable to the transverse component only in the strongly 
degenerate regime (see Figure 4), while the axial and mixed components 
(the latter contributes only to the difi^erential rates and emissivities but not 
to the total, since it is anti-symmteric in Ei and E2) are negligibly small. 
Notice that the first few Legendre coefficients are all comparable in mag- 
nitude. Moreover, the emission process is anisotropic (sec also Figure 7 of 
Ref. 16). "Light" photon (wp <C T) decays result in neutrino pairs with 
small outgoing angles between them, whereas "massive" photons {u)p T) 
yield back-to-back neutrino emission. 

Figure 4 shows the individual contributions to the total emissivity 
from the transverse, longitudinal, and axial channels at T = lO^^K and 
T = lO^K, respectively. The curves show results from expressions derived 
by exploiting the Lenard identity. The symbols "x" and "+" show re- 
sults obtained by integrations of the differential emissivities. At all densities 
and temperatures, the contribution of the axial channel to the total emis- 
sivity is negligible. For each temperature, the emissivity in the transverse 
channel dominates over that of the longitudinal channel at low densities. 
However, the peak values in these two channels are attained at nearly the 
same density; thereafter their individual contributions coincide. 

In terms of the density and temperature dependencies of the plasma fre- 
quency Up and the electron chemical potential /Zg (see Figure 4 of Ref. [16]), 
a detailed qualitative and quantitative analyses of the basic features of Fig- 
ure 4 are provided in Ref. [16]. Consider the case of Qt first. At a fixed 
temperature, the basic features to note are: 




xnB{uJ,T)5\K ~Qi-Q2){2T:f{\Mf) 




(17) 




(18) 
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Figure 3. The transverse part of the production kernels for T = 8.62 MeV and psYg < 
10^^ g cm"''. The Legendre coefficients are shown for 1 = through 3. The neutrino 
energies Ei and E2 are in MeV and the Legendre coefficients are in units of 10^^^ 
erg~® cm^ s"'^. 



(1) Qt is independent of the density psYg until a turn-on density (/OsVe)to 
is reached, 

(2) For densities larger than this turn-on density, Qt exhibits a power-law 

rise until a maximum is reached at {pBYe)pca.k, and 

(3) For psYe 3> {pBYe)pca.k, tlic fall-off with density is exponential. 

In the case of T 3> (non-degenerate plasma), the main contribution 
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Figure 4. Individual contributions from the transverse, longitudinal, and axial channels 
to the neutrino emissivity. The mass density of protons in the plasma, psYe = TUpUe, 
where rup is the proton mass. Ye = rie/ns is the net electron fraction (riB is the baryon 
number density) , and Tie is net electron number density. 



to Qt comes from high photon momenta. To a very good approximation,!® 



where C(3) ~ 1.202 is Riemann's Zeta function and tut is the transverse 
photon mass. 

For T <C (degenerate plasma) , Qt takes the form^® 



^ E/ (Cy) i5/2y3/2 -c.,/r 

ASn^a 



(20) 



The analyses of the longitudinal and axial cmissivities can be carried out 
along the same lines as that for the transverse emissivity. 
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3.2. The Photoneutrino Process 

The leading order diagrams for the photoproduction of neutrino pairs, + 
7* +fe,^,T + j^e,/i,r) are shown in Figure 5. The kernels of the collision 

integral for the photoneutrino process can be calculated along the same 
lines as for the plasma process and are detailed in Ref. 17. 



e t'e e e' 




7* e" 7* e" 



Figure 5. Leading order Feynman diagrams describing the emission of a neutrino pair 
from the photoneutrino process. The charged current W - exchange channel produces 

only fei'e, whereas the neutral Z - exchange results in pairs of all three neutrino 
(e, fj,, and r) flavors. Contributions from positrons are obtained by the replacement 
e~ ^ e+. 



Numerical results for the symmetric and anti-symmetric components of 
mEg,Eg,) and '^{{Eg,Eg,) are shown in Figure 6 for the case T = lO" 
K= 8.62 MeV and psYe = 1 g cm~^. The results explicitly show the ex- 
pected symmetry properties in neutrino energies Eg and Egi . A comparison 
of the relative magnitudes (^\{Eg,Eqi) (not shown here, but see Figure 6 
of Ref. [17]) and <^l''^^'^\Eg,Eg,) shows that the Z = term is the domi- 
nant term. The magnitude of $o f^^^™) amounts to only 10% of the leading 
<|)P(sym) gQj^|;j.j]3^|;jQj^ ^]-^(, contributions of i^^^^^^™-^ and ^p'^^^™) are 6% 
and 3%, respectively. In physical terms, this means that neutrino-pair emis- 
sion from the photo-neutrino process is dominantly isotropic. Therefore, 
depending on the required accuracy, $q might be adequate in practi- 
cal applications. Note also that that the production kernels are negligible 
for energies Eg and Egi > lOT. 
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Figure 6. Symmetric and anti-symmetric parts of the Legendre coefficients <l?g and 



f.P (asym) 



in the production kernels for the transverse case. The f-energies Eq and E„ 



are in MeV and the Legendre coefficients are in units of 10^^^ Tfi erg 



Contributions of the transverse and longitudinal components to the total 
emissivity are shown in Figure 7. For all temperatures at sufficiently high 
net electron densities rig (i.e., a degenerate plasma), the inequalities 

/ie > T, A^e > "T-e, 

UJP » T, Wp > TOe (21) 

are satisfied. In this case. 



Qi 



4 aGliCl + Cj 



(22) 



3 (27r)6 

The nondegenerate situation occurs at sufficiently low densities for 
which /ie — nie ^ T. In this case, both Qt and Ql exhibit a plateau 
for temperatures T > 10^ K. The emissivity and rate can be expressed in 
terms of the simple expressions 

^Qt\ 20aGl{Cl + Cl) /T^ x 775. 54^ 

Its X 136.50; 



(23) 



TtJ 3(27r)6 
The significance of the detailed differential rates and emissivities of neu- 
trino emission and absorption processes in calculations of core-collapse su- 
pernovae in which neutrino transport is strongly coupled with hydrody- 
namics remains to be explored. 



4. Protoneutron Stars 

A protoneutron star (PNS) is born in the aftermath of the gravitational 
collapse of the core of a massive star accompanying a successful super- 
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Figure 7. Individual contributions from the transverse and longitudinal channels to the 
neutrino emissivity as a function of baryon density at the indicated temperatures. The 
error bars show the variance of the Monte Carlo integration. For densities in excess 
of nuclear density shown by the dotted vertical line, neutrino production from strongly 
interacting particles dominate over QED-plasma processes. 



nova explosion. During the first tens of seconds of evolution, nearly all ('^ 
99%) of the remnant's binding energy is radiated away in neutrinos of all 
flavors. ^^'^^'^'^'^'''^^'^^ The neutrino luminosities and the emission timescale 
are controlled by several factors, such as the total mass of the PNS and 
the opacity at supranuclear density, which depends on the composition 
and dense matter equation of state (EOS). One of the chief objectives in 
modeling PNS's is to infer their internal compositions from neutrino sig- 
nals detected from future supernovae by SuperK, SNO and others under 
consideration, including UNO.^^ 
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4.1. The Evolution of a Protoneutron Star 

The evolution of a PNS proceeds through several distinct stages^^'^° and 
with various outcomes, as shown schematically in Figure 8. Immediately 
following core bounce and the passage of a shock through the outer PNS's 
mantle, the star contains an unshocked, low entropy core of mass ~ 0.7 
M0 in which neutrinos arc trapped (stage 1 in the figure). The core is 
surrounded by a low density, high entropy (5 < s < 10) mantle that is both 
accreting matter from the outer iron core falling through the shock and 
also rapidly losing energy due to electron captures and thermal neutrino 
emission. The mantle extends up to the shock, which is temporarily stalled 
about 200 km from the center prior to an eventual explosion. 



-20 hn 




idea t ill niri^t Ion) 
block hale 



cccrctiDn 
shock IHt-o 




(1) t = s 



deleplonization^ 
core healing 
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^jbservablc cs X-ray 
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Figure 8. The main stages of evolution of a protoneutron star. Shading approximately 
indicates relative temperatures. 



After a few seconds (stage 2), accretion becomes less important if the 
supernova is successful and the shock has ejected the stellar envelope. Ex- 
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tensive neutrino losses and deleptonization will have led to a loss of Icpton 
pressure and the collapse of the mantle. If enough accretion has occurred, 
however, the star's mass could increase beyond the maximum mass capable 
of being supported by the hot, lepton-rich matter. If this occurs, the rem- 
nant collapses to form a black hole and its neutrino emission is believed to 
quickly cease. 

Neutrino diffusion deleptonizes the core on time scales of 10-15 s (stage 
3). Diffusion time scales are proportional to i?^(cAy)~^, where R is the 
star's radius and is the effective neutrino mean free path. This generic 
relation illustrates how both the EOS and the composition influence evolu- 
tionary time scales. The diffusion of high-energy (200-300 MeV) u^s from 
the core to the surface where they escape as low-energy (10-20 MeV) i/'s 
generates heat (a process akin to joule heating). The core's entropy ap- 
proximately doubles, producing temperatures in the range of 30-60 MeV 
during this time, even as neutrinos continue to be prodiguously emitted 
from the star's effective surface, or i^— sphere. 

Strange matter, in the form of hyperons, a Bose condensate, or quark 
matter, suppressed when neutrinos are trapped, could appear at the end of 
the deleptonization. Its appearance would lead to a decrease in the maxi- 
mum mass that matter is capable of supporting, implying metastability of 
the neutron star and another chance for black hole formation.^^ This would 
occur if the PNS's mass, which must be less than the maximum mass of 
hot, lepton-rich matter (or else a black hole would already have formed), 
is greater than the maximum mass of hot, lepton-poor matter. However, if 
strangeness does not appear, the maximum mass instead increases during 
deleptonization and the appearance of a black hole would be unlikely unless 
accretion in this stage remains significant. 

The PNS is now lepton-poor, but it is still hot. While the star has 
zero net neutrino number, thermally produced neutrino pairs of all flavors 
dominate the emission. The average neutrino energy slowly decreases, and 
the neutrino mean free path increases. After approximately 50 seconds 
(stage 4), A ~ ii, and the star finally becomes transparent to neutrinos. 
Since the threshold density for the appearance of strange matter decreases 
with decreasing temperature, a delayed collapse to a black hole is still 
possible during this epoch. 

Following the onset of neutrino transparency, the core continues to cool 
by neutrino emission, but the star's crust remains warm and cools less 
quickly. The crust is an insulating blanket which prevents the star from 
coming to complete thermal equilibrium and keeps the surface relatively 
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warm (T w 3 x 10^ K) for up to 100 years (stage 5). The temperature 
of the surface after the interior of the star becomes isothermal (stage 6) 
is determined by the rate of neutrino emission in the star's core and the 
composition of the surface. 

4.2. Neutrino Signals in Terrestrial Detectors 




time (s) 

Figure 9. The evolution of the total neutrino luminosity for npQ PNS's. Shaded bands 
illustrate the limiting luminosities corresponding to a count rate of 0.2 Hz, assuming a 
supernova distance of 50 kpc for 1MB and Kamioka, and 8.5 kpc for SNO and SuperK. 
The widths of the shaded regions represent uncertainties in the average neutrino energy 
from the use of a diffusion scheme for neutrino transport. 

A comparison of the signals observable with different detectors is shown 
in Figure 9, which displays L^, as a function of baryon mass Mb for stars 
containing quarks in their cores. In the absence of accretion, Mb remains 
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constant during the evolution, while the gravitational mass Mq decreases. 
The two upper shaded bands correspond to estimated SN 1987A (50 kpc dis- 
tance) detection limits with KII and 1MB, and the lower bands correspond 
to estimated detection limits in SNO, SuperK, and UNO, for a Galactic 
supernova (8.5 kpc distance). The detection limits have been set to a count 
rate dN/dt = 0.2 Hz.^^ It is possible that this limit is too conservative 
and could be lowered with identifiable backgrounds and knowledge of the 
direction of the signal. The width of the bands represents the uncertainty 
in < Ej}^ > due to the diffusion approximation. ^^'^^'^^ It is possible to dis- 
tinguish between stable and metastable stars, since the luminosities when 
metastability is reached are always above conservative detection limits. 

4.3. Metastable Protoneutron Stars 

Protoncutron stars in which strangeness appears following dclcptonization 
can be metastable if their masses are large enough. One interesting diagnos- 
tic that could shed light on the internal composition of neutron stars would 
be the abrupt cessation of the neutrino signal. This would be in contrast 
to a normal star of similar mass for which the signal continues to fall until 
it is obscured by the background. In Figure 10 the lifetimes for stars con- 
taining hyperons (npH), kaons (npK) and quarks (npQ) are compared. 
In all cases, the larger the mass, the shorter the lifetime. For the kaon and 
quark PNSs, however, the collapse is delayed until the final stage of the 
Kelvin-Helmholtz epoch, while this is not necessarily the case for hyperon- 
rich stars. In addition, there is a much stronger mass dependence of the 
lifetimes for the hyperon case. 

Clearly, the observation of a single case of metastability, and the deter- 
mination of the metastability time alone, will not necessarily permit one to 
distinguish among the various possibilities. Only if the metastability time is 
less than 10-15 s, could one decide on this basis that the star's composition 
was that of npH matter. However, as in the case of SN 1987A, independent 
estimates of Mb might be available.^® In addition, the observation of two 
or more metastable neutron stars might permit one to differentiate among 
these models. 
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Figure 10. Lifetimes of metastable stars versus the PNS baryon mass Mb- Thick lines 
denote cases in which the maximum gravitational masses of cold, catalyzed stars arc 
near 1.45 Mq, which minimizes the metastability lifetimes. The thin lines for the npQ 
and npH cases are for EOSs with larger maximum gravitational masses (1.85 and 1.55 
M0, respectively.) 

5. Neutron Stars in Their Old Age 

See the contribution from Dany Page to this symposium. 

6. Outlook 

The advent of new-generation neutrino detectors such as Super- 
Kamiokande and the Sudbury Neutrino Observatory promises thousands of 
neutrino events in the next Galactic supernova. These will provide crucial 
diagnostics for the supernova mechanism, important limits on the released 
binding energy and the remnant mass, and critical clues concerning the 
composition of high density matter. Research in this area will ascertain 
the extent to which neutrino transport is instrumental in making a super- 
nova explode. Other bonuses include the elucidation of the possible role of 
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supernovae and neutrinos in r— process nucleosynthesis. 

The main issues that emerge from PNS studies concern the metasta- 
bihty and subsequent coUapse to a black hole of a PNS containing quark 
matter, or other types of matter including hyperons or a Bose condensate, 
which could be observable in the v signal. However, discriminating among 
various compositions may require more than one such observation. This 
highlights the need for breakthroughs in lattice simulations of QCD at finite 
baryon density in order to unambiguously determine the EOS of high den- 
sity matter. In the meantime, intriguing possible extensions of supernova 
and PNS simulations with npQ and npK matter include the consideration 
of heterogenoeus structures and quark matter superfluidity. See also the 
contribution from Sanjay Reddy to this symposium. 

7. Conclusions 
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